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SUMMARY 


A reduced-order control law is synthesized by minimizing a performance 
index defined by a weighted sum of mean-square steady-state responses and con- 
trol inputs, as in a stationary linear quadratic Gaussian (LQG) analysis. The 
order of the control law is assiomed to be less than the order of the plant. 
Gradients of the performance index with respect to the design variables in the 
reduced-order control law are determined by solving a pair of Lyapunov equa- 
tions. An analogy with the LQG solution is utilized to select a set of design 
variables and their initial values. Using the gradients, a nonlinear program- 
ing algorithm searches for the control law design variables which minimize the 
performance index. During the design cycle, an input-noise adjustment pro- 
cedure is introduced to improve the system's stability margins. Thus an 
optimal, reduced-order, robust feedback control law can be synthesized without 
truncating the plant model order. 

This control law synthesis method was applied to the synthesis of an 
active flutter- suppress ion control law for a wind-tunnel model of an aero- 
elastic wing represented by a 25th-order plant. The resulting fourth-order 
control law is shown by analysis to increase the flutter dynamic pressure by 
at least 44 percent with good stability margins, while minimizing the control 
input. Although theoretically proven only for a full-order controller, the 
input-noise adjustment procedure also improved the phase and gain margins for 
the reduced-order controller. Fourth-order control laws were also obtained 
by truncation and residualization methods, and results of these control laws 
were compared with those of the present reoptimization method. The study 
indicates that by using the present algorithm, nearly optimal low-order con- 
trol laws with good stability margins can be synthesized. 

INTRODUCTION 

The state-space equations describing control problems involving flexible 
structures are usually of high order. This characteristic is particularly true 
for an aeroelastic system (plant) which requires a large number of states to 
accurately represent the flexible structure, unsteady aerodynamics, and actuator 



dynamics (refs. 1 to 3) . The order of a realistic design problem could be 60 or 
more (refs. 1 and 4) . An optimal feedback control law based on the standard 
linear quadratic Gaussian (LQG) solution would be of the same high order as the 
plant. This control law (controller) is usually very sensitive to modeling 
errors, has poor stability margins, and often is too complex to implement in a 
flight computer. In this paper, a method is developed using optimization tech- 
niques for designing a reduced-order control law without these disadvantages. 

Reduced-order control laws have been designed using transfer function 
matching, modal truncation, and residualization techniques (refs- 4 to 6) . 

These methods only approximate the full-order controller, and the resulting con- 
trol law is no longer optimal. In the present approach, the control law is 
synthesized by minimizing a performance index defined by a weighted sum of mean- 
square steady-state responses and control inputs as in a stationary LQG analysis. 
However, the order of the control law is assumed to be less than the order of 
the plant. Gradients of the performance index with respect to the design vari- 
ables of the reduced-order control law are determined by solving a pair of 
Lyapunov equations. Using the gradients, a nonlinear programing algorithm 
searches for the design variables which minimize the performance index. The 
basis of this reduced-order controller design procedure is described in refer- 
ences 7 and 8 and has been applied in reference 9 for attitude control of a 
flexible spacecraft. 

Two problems arise when applying this general method: selection of a set 
of design variables and their initial values which result in a stable system. 

In the present paper, a new systematic approach is developed to overcome both 
problems. The design variables of the reduced-order controller are chosen so 
that they are analogous to the optimal full-state feedback and Kalman estimator 
gain matrices. These full-order optimal gain matrices are used for choosing the 
initial values of the design variables. In the limit when the order of the con- 
troller is the same as the order of the plant, the algorithm provides a solution 
which is identical to the optimal LQG results. 

In general, plants using observer-based control laws exhibit poor stability 
margins. Recently, Doyle and Stein (ref. 10) presented an input-noise adjust- 
ment procedure to improve the stability margins of full-order LQG control laws. 
This procedure is used in the present design algorithm. Results indicate that 
the procedure can improve the stability margins of the system with a reduced- 
order control law, although the theoretical proof exists only for a full-order 
control law. 

The method of this paper is then applied to the synthesis of a 4th-order 
flutter- suppression control law for a wind-tunnel model of an aeroelastic wing, 
represented by a 25th-order plant. The performance of the flutter-suppression 
system (FSS) using the reduced-order (4th-order) control law is analyzed and 
compared with the perfoirmance using the corresponding full-order (25th-order) 
control law. A comparison of results using fourth-order control laws obtained 
by truncation and residualization methods is also presented. 
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EQUATIONS OF MOTION IN STATE-SPACE FORM 
Plant Model 

The equations of motion for an aeroelastic system can be written (refs. 2 
and 3) as 


([m]s2 + [dJs + [k]){?p) + q[Q 


r-1 


r-“l r?r 


GIV 


( 1 ) 


where 


[m] 

[k] 


generalized mass matrix 


structural damping matrix 


generalized stiffness matrix 


Laplace variable 




generalized coordinate vectors for flexible modes and 
for control surface deflections 


gust velocity vector 
dynamic pressure 
free-stream velocity 


(The symbols used in this report are listed after the references on page 36.) 


The s-plane approximations of the unsteady aerodynamic forces 
are expressed as 


and 


= H - [y(n) * H(f)' 




m=l ( s + 3 

' c m 


( 2 ) 


where 


reference chord length 


3 



L 


nuinber of aerodynamic lag terms 


3m aerodynamic lag 


The real coefficient matrices 




relate unsteady aerodynamics 


the frequency plane to the approximating functions in the s-plane. Unsteady 
aerodynamic forces on the wing and control surfaces due to sinusoidal motion 
and gust are generated herein using doublet lattice aerodynamics (ref. 11) . 
denoting the flexible modes, their time derivatives, and the aerodynamic lag 
terms within the summation sign in equation (2) by the state variable vector 
equation (1) and an accelerometer sensor output equation can be expressed in 
state-space form (see appendix A) as 


in 


By 

X*, 


{x*} = [f*]{x*} + I G 2 j{u*} 

(3a) 

{y} = [h^J{x*} + |^H 2 ' H 3 j{u*} 

(3b) 


where 


F* 

* * 
^1 ^ ^2 


dynamics matrix (see eqs. (A7) ) 
input matrices (see eqs. (A7) ) 
measurement matrices (see eqs. (A9) ) 


and the state vector and input vector are 



where x* are aerodynamic states (eq. (A3)). As shown in appendix A, the 
^m 

derivative terms in u* can be eliminated by augmenting x* with a model of 
the control surface actuator dynamics and a model of the turbulence. The state- 
space model of the plant can then be written in the standard form as 


Xs = FXg + G^u + G„w 


y = Hx_ + V 


= Vs 


(4) 


(5) 

( 6 ) 
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where 


Xg augmented plant state vector composed of x*, actuator states, and 

turbulence states (Ng x l) 

u control input vector (Nq x i) 

w plant noise vector (N^ x i) 

y measurement output vector (Nq ^ 1) 

V measurement noise vector (Nq ^ D 

design output vector (Nq x i) 


and the matrices in equations (4) to (6) are defined in equations (A16) , (A17) , 
and (A19) . The noise vectors w and v are modeled as zero-mean white noise 
processes with intensity matrices and respectively. 

Controller Model 

A block diagram of the control scheme is shown in figure 1. The controller 
model is assumed to be 

x^ = AXc + By (7) 

u = CXc (8) 

where x^ represents the controller state vector of order M where M < Ng. 

The terms A, B, and C are controller dynamics, input, and output matrices of 
size M X M, MX and N^ ^ M, respectively. Of the M(M + Nq + N^) ele- 

ments of the A, B, and C matrices, only M(Nq + N^) are independent 
(ref. 12) and can be chosen as free design variables- Equation (7) represents 
a filter which processes the output measurements before being fed back through 
the gain matrix C. Together, equations (7) and (8) represent a transfer func- 
tion relation u = c\sl - a] By which is commonly referred to as an output 
feedback control law. 


CONTROL LAW SYNTHESIS METHOD 

The control law is synthesized by minimizing a performance index defined by 
a weighted sum of mean-square steady-state responses and control inputs, as in a 
stationary LQG analysis. However, the order of the control law is assumed to be 
less than the order of the plant. Gradients of the performance index with 
respect to the design variables in the reduced-order control law are determined 
by solving a pair of Lyapunov equations. An analogy with the LQG solution is 
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utilized to select a set of design variables and their initial values. Using 
the gradients, a nonlinear programing algorithm searches for the control law 
design variables which minimize the performance index. During the design cycle, 
an input-noise adjustment procedure is introduced to improve the system's sta- 
bility margins (robustness) . 


Augmented State Equations 


By defining an augmented state vector 



the closed-loop system is represented by 


or 



X 


a 




+ G^n 


( 9 ) 


( 10 ) 


If Fq is stable with the chosen values of A, B, and C, the covariance of 
Xg reaches a steady-state value which satisfies the Lyapunov equation given by 




( 11 ) 


where 




lim El 

t-K» 


r - 

1 

X 1 

^SC 

^a^a = 

„ T ' 

X 

L J 

^SC 1 

c 


( 12 ) 


and 


Ra — 


R ' ^ 

_w I 0_ 

0 , ^ 


( 13 ) 
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Performance Index 


The control law is synthesized using a conjugate gradient algorithm 
(ref. 13) to search for the M(Nq + N^) design variables of the Mth-order 
control law (eqs. (7) and (8)) which minimize a performance index J defined 
by 



where and Q« are symmetric weighting matrices and X and U are 

JL ^ s 

steady-state covariance matrices of the states Xg and control inputs u, 

respectively- Substituting equation (8) into equation (14) and using the 

definitions given by equations (9) and (12) result in 


J 


tr 




__q__ 

T 

C Q 2 C 


sc 




(15) 


Gradients 

It is shown in appendix B that the gradients of J with respect to the 
elements of the matrices A, B, and C can be expressed as 


9j 

3a 


2 



+ A^x 


c 


(16) 


8j 

8b 


2 



+ A^X 






(17) 


9j 

8C 




2 Q2CXc 


+ G, 




■^s^sc 


^sc^c 


)] 


(18) 
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L 


where A is a (Ng + M) x (Ns + M) symmetric Lagrange multiplier matrix 
defined as 


A E 




S 1 

“SC 

iFi 

1 

1 

^c 

SC j 



( 19 ) 


which satisfies the dual Lyapunov equation 
fJa + AFa = -Qa 


( 20 ) 


Design Variables 

There are several possible ways of selecting a set of M(Nq + N^) free 
design variables of the control law (eqs. (7) and (8)). In a canonical form 
of the control law, the coefficients of numerator and common denominator poly- 
nomials can be identified as the design variable set. In a block diagonal 
form, the real and imaginary parts of the controller poles and numerator resi- 
dues can be taken as the design variable set. While these representations can 
be incorporated into the present algorithm (to optimize an existing control law 
in transfer function form) , they are not very general and are often inconvenient 
from a matrix bookkeeping point of view. Also the initial values of the design 
variables can be difficult to estimate- In the present method the M(Nq + N^) 
elements of B and C are chosen as the free design variables and a subset of 
the Kalman estimator gain and optimal full-state feedback gain matrices obtained 
from the LQG solution are chosen as their initial values. The justification is 
as follows . 

If for a full-order controller (M = Ng) , we let A = F - BH + G^C and use 
the elements of B and C as design variables, then it is shown in appendix C 
that the optimized values of B and C are identical to the steady-state 
Kalman estimator gain and optimal full-state feedback gain matrices, respec- 
tively. The controller states x^ are identified as the asymptotic estimates 
of the plant states Xg. By analogy, when M < Ng, the low-order controller 
may be treated as a partial estimator of M key states denoted by x^^ E RXg 
(R is a Boolian matrix of order M x . a full discussion of this approach 
and the asymptotic behavior of the estimation error is presented in appendix D. 

Thus we initially set 


A = 

R(F - BqH + G^Cq)r'^ 

(21) 

B = 

RBq 

(22) 

C = 

o 

o 

(23) 
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where M key states are chosen for estimation. The matrices Bq and Cq are 
the full-order Kalman estimator and optimal full-state feedback gain matrices, 
respectively. Siibsequently , 

A = - BHR^ + RGuC (24) 


Since A is a function of B and C, the total gradients of J with respect 
to B and C are given by 


^ ^ T 

dB 9 b 3a ^ 




dJ _ dJ 
dC 9C 



T ^ 
9a 

J 


(25) 


Svibstituting the partial gradients from equations (16) to (18) in equations (25) 
results in 




dJ 

dB 


= 2 


■^SC ^S ■^C^SC (^sc ^sc ^ ^c^c ^ -^C^^ 


dJ 

dC 


= 2 


+ <^u{^s^sc + ^sc^c) + 


SC c c 


(26) 


With the performance index J and its gradients known, a conjugate gradient 
optimization procedure (ref. 13) can be used to minimize J. The (Ns + M) th 
order Lyapunov equations are solved by a program described in reference 14. 
This process is feasible only if is a stable matrix for each iteration. 


Robustness Recovery Technique 


Plants using full-order LQG controllers usually exhibit poor stability 
margins. It has been shown in reference 10 that by introducing a fictitious 
input noise r|^ of intensity R^ in equation (8) (the noise is not included in 
the performance index directly) and by increasing Ru by a scalar factor, the 
full-order controller system loop transfer function (loop broken at u) asymptot- 
ically approaches the corresponding full-state feedback loop transfer function 
(and hence they have the same stability properties) . This input-noise adjust- 
ment procedure is incorporated in the present reduced-order controller design 


process by replacing 
,T 

, respectively. 


and 


T I m 
Bu I ^ 


and w in equation (10) by pu I 

Consequently, R^ in equation (13) is replaced by 
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I 0 

_0 


. . T 

The resulting effect is an additional noise teim G^Ru^u first Ng x Ng 

sets of Lyapunov equation (11) . In appendix D, it is shown that r)u appears as 
an additional noise parameter in the partial state estimation error dynamics. 
Although not proven theoretically/ numerical results show that this input-noise 
adjustment procedure is also capable of improving stability margins of the plant 
with a reduced-order controller. The input-noise intensity matrix becomes 

a major design parameter in this controller synthesis technique. 


Design Algorithm 

A block diagram of the design algorithm is shown in figure 2. The design 
algorithm consists of the following steps: 


1 . 


2 . 

3. 


4. 


5. 


6 . 
7. 


8 . 


9 . 


Obtain 

Bq and Cq . 


Select 

M key states (matrix 

R) . 

Compute 

T T 

RFR , RG^, HR , B 

T 

= RBq, and C = CqR . 

Set A = 

= RFR*^ - BHR*^ + RG^C. 


Test Fg 

for stability. 


If Fa 

is stable, solve Lyapunov equations (11) and 


and A. 

Compute J, dJ/dB, and dJ/dC. 

Obtain the next set of B and C by a conjugate gradient algorithm. 

Repeat steps 4 through 8 with new values of B and C until J 
converges to a minimum value. 


In step 2, it is often convenient to transform the plant state equations to a 
block diagonal form for easier selection of key states to be estimated. If 
is initially unstable in step 5, then either select a new set of states or 
select a new value of R^ and return to step 1/ until a stable is obtained. 

If F^ becomes unstable during a linear search in the optimization process, the 
algorithm automatically searches in a new direction from the last stable F^ 
solution. 


APPLICATION TO FLUTTER SUPPRESSION 
Model Description 

The control law synthesis method described in the previous sections is 
applied to the synthesis of an active flutter-suppression control law for an 
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aeroelastic wind-tunnel model. The geometry of the cantilever wing model, along 
with the control surface and sensor (accelerometer) locations, is shown in fig- 
ure 3. The model is described in detail in reference 1. The state-space equa- 
tions are derived using five flexible modes and two aerodynamic lag terms (L = 2, 
3i = 0.225, 32 “ 0.5). All structural damping is assumed to be zero. The 
actuator dynamics are represented by a third-order transfer function given by 
(ref. 1) : 




C 




(214) (89 450) deg 

(s + 214) (s2 + 179.45s + 89 450) 


(27) 


A turbulence model represented by the following second-order transfer function 
is used to approximate a Dryden gust spectrum: 



a 

wg 


root-mean- square (rms) gust velocity, m/sec 


Z scale of turbulence or characteristic length, m 

V flight velocity, m/sec 


(28) 


For numerical calculation £ = 30.48 m and O =1. The resulting plant model 

wg 

is a single-input single-output system of order 25. The input u is the actu- 
ator command signal 6^ and the output y is the accelerometer signal z. The 
dynamic-pressure root locus at Mach 0.9 is shown in figure 4 without the flutter- 
suppression system. Flutter is predicted at a dynamic pressure of 4.7 kPa and a 
frequency of 50 rad/sec. The present objective is to design a low-order control 
law which will increase the maximum operating dynamic pressure by 44 percent 
above the experimental flutter dynamic pressure of 5.32 kPa (ref. 1), with mini- 
mum control surface activity and adequate stability margins. The control law 
design point is thus chosen to be a dynamic pressure of 7.66 kPa at Mach 0.9. 


Control Laws 

Ful l-order control law .- For comparison purposes, the full-order control 
law, which is identical to the LQG solution, is obtained first. The design out- 
put matrix is chosen to be the same as the sensor output matrix H. The 

weighting matrices, which are scalars for this single-input single-output sys- 
tem, are selected to be Qi = 0.0001, Q 2 = 50 000, = 1.0, and = 1.0. 
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This choice of weighting nearly reflects the unstable pole about the imaginary 
axis without affecting the stable poles. The unit plant noise intensity matrix 
results in unit rms gust velocity in equation (28) . The control laws are 
designed at three values of Rq/ namely 0, 0.00001, and 0.0001. Numerical val- 
ues of all these parameters were chosen on the basis of some preliminary results 
obtained with a 20th-order plant model (ref. 15) . To determine the improvements 
in the control law designed with increasing the phase and gain margins of 

the plant plus the control law are determined. The Nyquist diagrams of the 
plant with the full-order output feedback control* law, presented in figure 5, 
show a progressive increase in phase and gain margins. As predicted in refer- 
ence 10, the Nyquist diagram asymptotically approaches the full-state feedback 
result which is a circle of unit radius centered at (-1,0). The price paid is 
an increase in the optimal design performance index J and associated increases 
in rms values of the input and output (table I) . The closed-loop dynamic- 
pressure root locus using the full-order control law designed with Rq = 0.00001 
is presented in figure 6. The flutter dynamic pressure is increased to 9.0 kPa. 
Beyond this point, the control law complex pole (filter mode) is unstable. 

Fourth-order control law .- Fourth-order control laws were synthesized by 
selecting the first two flexible modes and their time derivatives as the four 
key states. The initial values of the gain matrices B and C were obtained 
from the full-order results. In transfer function form, the control laws 
designed with the three values of R^ are, respectively 

u ^ (-364.4) (s - 136.4) (s^ + 73.69s + 5697) 

z (s + 2.057) (s + 2057) (s^ + 46.37s + 2047) 

u ^ (1939.4) (s + 24.74) (s^ + 87.63s + 13 806) 
z (s + 3.864) (s + 3270) (s2 + 20.97s + 1423) 

u ^ (1424.6) (s + 31.90) (s^ + 100.0s + 12 210) 
z (s + 1.798) (s + 2541) (s2 + 17.76s + 1547) 


deg/g unit (r^ = o) (29) 


deg/g unit (r^ = o. 00001) (30) 


deg/g unit (R^ = O.OOOl) (31) 


Figure 7 shows a typical convergence pattern of the performance index 
obtained in the course of synthesizing control law (eq. (30)). Convergence is 
assumed when the change in the performance index is less than 0.1 percent of its 
initial value for three successive iterations. Most of the reduction in the 
performance index is achieved in the first few iterations. The optimized per- 
fomance index is within 16 percent of the full-order controller performance 
index . 

The Nyquist diagrams of the plant plus fourth-order controllers given by 
equations (29) to (31) are presented in figure 8. With the control law designed 
without input noise (eq. (29)), the stability margins are poor and relatively 
unaffected by controller order reduction. Although theoretically proven only 
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for a full -order controller, it is interesting to note that the input-noise 
adjustment procedure of reference 10 is also capable of improving the phase 
and gain margins for a reduced- order controller. For a given reducing the 

controller from full order to fourth order results in only a slight degradation 
(increase) in the optimal design performance index J (table I) . 

As a compromise between better stability margins and lower control surface 
motion (deflection and deflection rate) , the controller given by equation (30) , 
i.e./ designed with = 0.00001, was selected for overall performance analysis 

and comparison- The dynamic-pressure root locus using the fourth-order optimal 
control law of equation (30) is shown in figure 9. The closed- loop flutter 
dynamic pressure is 8.7 kPa. Beyond this point, the controller complex pole 
(filter mode) is unstable. A Bode plot of the plant plus this control law is 
presented in figure 10 and compared with the corresponding plot for the full- 
order case. The plots compare quite well up to approximately 400 rad/sec. A 
sharp drop in amplitude at 400 rad/sec for the fourth- order case is due to a 
zero in the plant (z/6c) transfer function. For the full-order case, this sharp 
drop in amplitude is removed by cancellation of the plant zero by a controller 
pole. 

The lowest order of the controller obtained in the course of this study is 
four. Higher order controllers are possible but do not result in substantial 
performance improvement. 


CONTROL LAW PERFORMANCE 

In this section, the performance of the fourth-order controller designed 
with R^ = 0.00001 (eq. (30)) is discussed and compared with the corresponding 
full-order (LQG) controller. 


Flutter Characteristics 

The dynamic-pressure root locus diagrams of the plant with the full 25th- 
order controller and 4th-order controller (eq. (30) ) are presented in figures 6 
and 9, respectively. The first two modes and the controller complex pole 
(filter mode) are less damped for the fourth-order controller than for the full- 
order controller, at dynamic pressures of 7.66 kPa and lower. The higher modes 
remain relatively unaffected by both control laws. With the full-order control- 
ler, the flutter dynamic pressure is increased by 91 percent (from 4.7 kPa to 
9.0 kPa) , whereas with the fourth-order controller, the increase is 85 percent 
(to 8.7 kPa) . Thus, the flutter-suppression performance deterioration due to 
this reduction in controller order is only 6 percent. In both cases, beyond 
the flutter point, the controller complex pole is unstable at a frequency of 
40 rad/sec . 


Stability Margins 

To examine the stability margins of the plant plus the control laws at off- 
design dynamic pressures, the phase and gain margins of this single-input 
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single-output system were computed over a range of dynamic pressures up to the 
design value (7.66 kPa) and are presented in figure 11. With the fourth-order 
controller at the design dynamic pressure, the gain margins are -5 dB and 
+12.3 dB. The phase margins are -46° and +53°. The closed- loop system is 
stable over the entire range of dynamic pressures, and a ±5 dB gain margin and 
±30° phase margin are maintained over the dynamic -pres sure range from 5.7 kPa 
to 7.66 kPa. This provides a measure of insensitivity of the designed active- 
control system to dynamic-pressure variation. Except for the positive phase 
margins which fall below +30° at dynamic pressures of 5.7 kPa and lower 
(fig. 11 ), all the other margins improve with decreasing dynamic pressure. 

Since the stability margins with the full-order controller also exhibit almost 
the same characteristics (shown by dashed line in fig. 11 ) , the deterioration in 
positive phase margin is not due to the controller order reduction. 


Steady-State Response 

Table I presents the steady-state root-mean- square (rms) responses of the 

control input control surface deflection 6 j^g, deflection rate 6 ^ 3 , 

and sensor output Z 2 mis design dynamic pressure using the fourth-order 

controllers (eqs. (29) to (31)) and corresponding full-order controllers. These 
values were obtained by setting = 0 and executing the usual covariance 

analysis with the optimized control law. With the fourth-order controller 

designed with R^ = 0.00001 (eq. (30)), and are 17.0° and 

756.0 deg/sec, respectively (all calculations performed for a 1 m/sec rms gust 
velocity) . The respective values with the corresponding full-order controller 
are slightly higher. These values are in the vicinity of the typical maximum 
allowable limits. 

To examine the steady-state ^rms' ^rms' ^rms off-design dynamic 

pressures, their values normalized by their design-point values are presented 
in figure 12 for the fourth-order controller designed with R^j = 0.00001 
(eq. (30)). These responses are computed with R^ = 0 and R^ = 0 and are 
slightly lower than those with R^ = 1.0. Figure 12 indicates that for dynamic 

pressures below the design point, films' ^rms ' ^rms consistently 

lower than their design-point values except for z^^g at 6.7 kPa which is 

2.6 percent higher. The deterioration of positive phase margin at lower dynamic 
pressures (fig. 11 ) does not adversely affect the system steady-state response. 


Transient Response 

The present design technique is based on minimizing the system input and 
output at steady state. Therefore, it is important to check the transient 
response with the designed controllers. The state transition matrix technique 

(ref. 14) was used to determine 6 , 6 , and z responses to a 1 ° step pulse of 
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0.01-second duration as an input command for the full- and fourth-order control- 
lers at the design dynamic pressure. These responses are shown in figure 13. 

The responses using the fourth-order controller are quite similar to those 
using the full-order controller and decay rapidly in a few cycles . The peak 

— 

values are 6=1°, 6 = -200 deg/sec, and z = 0.6g, respectively- These are 

somewhat higher than values for the full-order controller. The lower damping of 
the plant with the fourth-order controller is consistent with the earlier obser- 
vation from the root-locus diagrams (figs. 6 and 9). 


Perfoirmance With 65th-Order Plant Model 

The 4th-order controller of equation (30) was designed using a 25th-order 
state-space model of the plant obtained using five flexible modes and two aero- 
dynamic lag terms- In order to examine the performance of this control law with 
a more accurate state-space model of the plant, it was applied to a 65th-order 
plant model obtained by including the first 10 flexible modes and 4 aerodynamic 
lag terms (L =4, = 0-2, 32 = 0.4, 3s = 0.6, 34 = 0.8). The actuator 

dynamics and turbulence models (eqs. (27) and (28)) remained the same. The 
Nyquist diagram of this 65th-order plant plus the 4th-order controller of equa- 
tion (30) at the design dynamic pressure is presented in figure 14. The gain 
margins are -5.4 dB at 62 rad/sec and +12.0 dB at 310 rad/sec . The phase mar- 
gins are -62^ at 45 rad/sec and +57° at 82 rad/sec. These margins are slightly 
better than those with the 25th-order plant model (fig. 8(b)). The stability 
margins at lower dynamic pressures (not presented) are similar to those of fig- 
ure 11- The dynamic-pressure root locus of the 65th-order plant plus the 4th- 
order controller is presented in figure 15. All the modes not included in the 
original 25th-order plant used for the controller design are found to be stable. 
The flutter dynamic pressure is 8.85 kPa. Beyond this point, the controller 
complex pole (filter mode) is unstable at a frequency of 40 rad/sec. 


COMPARISON WITH TRUNCATED AND RESIDUALIZED CONTROL LAWS 

Reduced-order control laws can also be obtained from optimal control theory 
results by truncation and residualization methods. In the truncation method, 
the low-order controller is obtained by partitioning the optimal full-state 
feedback and Kalman estimator gain matrices and retaining only the part of the 
optimal controller corresponding to the key states and their first derivatives. 
In the residualization method (ref. 6) , a first-order correction is added by 
retaining the static part of the controller corresponding to the rest of the 
states. The present method can be considered a reoptimization of the truncated 
control law. For a comparison of the performance of the control laws obtained 
by these three related methods, 4th-order control laws are also designed by 
truncation and residualization methods using identical base data and noise 
intensities (R^ = 0.00001) for the 25th-order plant at the design dynamic pres- 
sure. The Nyquist diagrams of the plant plus the full-order controller and of 
the plant plus 4th-order truncated, residualized , and reoptimized control laws 
are presented in figure 16. The corresponding stability margins are presented 
in table II. For the truncation method, the gain margins are -10.3 dB and 
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+6.5 dB, and the phase margins are -17^ and +57^. For the residualization 
method, the gain margins are -7.1 dB and +6.0 dB, and the phase margins are 
-55^ and +38°. The shapes of the corresponding Nyquist diagrams (figs. 16(b) 
and (c)) indicate that the phase and positive gain margins would decrease 
rapidly with increase in gain. In comparison, for the present method, the over- 
all stability margins are better and closer to the full-order optimal result. 
Since the corresponding Nyquist diagram is nearly circular around (-1,0) 

(fig. 16(d)), the stability margins are less sensitive to change in gain. 

The root-mean- square (rms) values of the steady-state responses films' 

^rms' ^rms using the control laws designed by the three methods are pre- 

sented in table II along with the corresponding values using the full-order 
optimal LQG controller, for comparison. These values are computed using Ru = 0 
and correspond to the steady-state response due to a 1 m/sec rms gust. The 

present method yields the lowest ^rms' while the lowest S^ms obtained by 

the residualization method. The truncation method provides the lowest 

In general, except in the truncation method, the control surface rms responses 
are of the same order. 


CONCLUSIONS 

A method of synthesizing reduced-order optimal feedback control laws for a 
high-order system is developed. A design algorithm which employs a nonlinear 
programing technique and an analogy with the LQG solution is presented. An 
input-noise adjustment procedure is used to improve the stability margins of 
the system. This general method is applied to the synthesis of an active 
flutter-suppression control law for a wind-tunnel model of an aeroelastic wing. 
The important results of this study are the following: 

1. The present method can be used to synthesize optimal, reduced-ord^r , 

robust feedback control laws for a high-order plant without truncating 
the plant model order. It is applicable to a multi-input, multi-output 
system. 

2. Application of this method to a 25th-order plant representing an aero- 

elastic wing model provides a 4th-order flutter-suppression control 
law which is shown by analysis to increase the flutter dynamic pres- 
sure by at least 44 percent (at Mach of 0.9) with good stability mar- 
gins while minimizing the control input. 

3. The numerical results indicate that the input-noise adjustment procedure 

of Doyle and Stein is capable of improving the stability margins of 
the system with a reduced- order controller, although the theoretical 
proof exists only for a full- order controller. 
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4. The present reoptimization algorithm appears to provide better control 
laws than those obtained by truncation and residualization methods. 


Langley Research Center 

National Aeronautics and Space Administration 
Hampton, VA 23665 
June 17, 1981 
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APPENDIX A 


STATE-SPACE EQUATIONS 

The detailed derivation of equations (3) to (6) are presented in this 
appendix - 


Plant Equation 


Let us define the vectors 



and also define the following matrices 



(Al) 


(A2) 


L) (A3) 


(A4) 


Separate the control surface and flexible mode components by partitioning the 


the matrices 


and 




as 


pn] - [[^n]^ j [ah]J 


(n = 0, 1, 2) (A5) 
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[%n] = [[^]f i [Bm]c] (m = 1, 2, . . L) (A6) 

Then by using the definitions in equations (Al) to (A6) , equation (1) can be 
written as shown in matrix equation (A7a) on page 21. With zero initial condi- 
tions assumed, equation (A7a) is denoted by equation (3a) in the time domain: 

{x*} = + [g* ; (A7b) 


Sensor Output Equation 


An accelerometer sensor output equation can be expressed as 

{v> - wQ 


(AS) 


where is a matrix of modal amplitudes at the sensor location- Using the 
definition in equation (A2) and the second row of equation (A7a) , equation (AS) 
can be written as 




(y) = M ^ [k] I [s] ; [i] ; [i] I > . . • w]< 




+ q 


m[«J 


-1 


■pole : I ; pi]j2v) ; p2]j2v)_ 




<- 


(A9a) 
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Matr ix equation ( A7a ) 


APPENDIX A 



APPENDIX A 


or 


{y} = |H*j {x*} + ^H* I 

Thus, equation (3b) is obtained. 


(A9b) 


Actuator and Turbulence State Equations 
The actuator transfer function can be expressed in state-space form as 


{*Jo> ' [''Sc] tic} * [®Sc](“) 




fQ 




> ■ fSc] ('‘Sc) 






(AlO) 


where A* , B* , and C* are the dynamics matrix, the control distribution 
ac ac ac 

matrix, and the output matrix of the actuator and {u} is the actuator command 
input. The order of the actuator state vector is the order of the denominator 
in the actuator transfer function. A model of the turbulence is used whose 
power spectral density approximates the power spectral density of the gust 
velocity. The model of the turbulence can be expressed as a Markov’s process 


ti} ■ * fa]'"’"] 


V J 







(All) 


where A*, B*, and C* are the dynamics, control distribution, and output 

matrices of the turbulence and {w} is a white noise process . 
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Plant State Equation 

Equations (AlO) and (All) can be combined and written as 



(A12) 



(A13) 


By augmenting the vector 
tion (A12) into equation (A7b) 




{x*} with 


ac 


and substituting equa- 


equations (3a) and (A13) can be written as 




F* 

1 

1 

G*C* 

1 ac 

1 

1 

2 G 

F = 

0 

1 

A* 

ac 

1 

0 


5 ” 

1 

0 

1 




1 


1 



(A14) 


(A15) 


(A16) 
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0^' 


oj 

G E 


G — 

5“ 

u 

ac 

w 



”5“ 


B* 




G 


one can denote equation (A14) by equation (4) . 

Sensor Output Equation in State Variables 
Substituting equation (A12) into equation (A9b) , one obtains 


{y} = 


I 


H* ' H*C* ‘ 
1 I 2 ac I 



+ {v} 


where a white noise process {v} is added to represent measurement 
uncertainties . 

Defining 




= IH* ' 

X I 


H*C* • H*C* 
2 ac , 3 G_ 


and using the definition of the state vector in equation (A15) 

in equation (5)- In a similar manner any design output vector 

to obtain equation (6). 


written in terms of 


(A17) 


(A18) 


(A19) 

results 
can be 
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GRADIENT OF PERFORMANCE INDEX 

Minimization of the performance index J defined by equation (15) in which 
the covariance matrix X^ satisfies the Lyapunov equation (11) can be treated 
as a constrained optimization problem. Equation (15) can be written as 


J = tr^Q^Xa + 


+ X^F, + G^R 
a a a a a a a 


n] 


(Bl) 


where A is a Lagrange multiplier matrix. Define a matrix P as 


P = 


0 I C 
B I A 


(B2) 


tl (Nc+M)x(No+M) 


/\ 

The necessary conditions for minimization of equation (Bl) with respect to P 
are then equation (11) (from 8J/9A = 0) along with 


0 = 


9j 

9X;= 


T T 

Q + F A + AF^ 
^a a a 


(B3) 


0 = 


9j 
9p. . 

ID 


= tr 


9Q 


^ X + a"" 
3. 


3p. . 

ID 


9f 


3f 


X^ + Xa 


a +_A 


9p. . 

ID 


9p. . 9p. . 

ID ID 


(caRaGa’’) 


(B4) 


where P^^ is a nonzero element of P. Since is symmetric, equation (B3) 

leads to Lyapunov equation (20) which is dual to equation (11) and can be solved 
to obtain A. Thus, A is also symmetric. Equation (B4) represents 

M(M + Nq + N^) equations for 9 j/9p. .. Since equation (B4) cannot be solved 
/\ ^ ^ 

explicitly for P, one can use it to express the gradients of J with respect 
/\ 

to P^j as follows- Since and A are symmetric, the trace proper- 

ties, of compatible matrix products can be used to write equation (B4) as 


9j 

/s 

9p. . 
ID 


= tr 


9q 9f 9g 

— X + 2 — X A + 2 — R G ^A 


9p. . 
ID 


9p. . 
iD 


9p. . 
ID 


a a 


(B5) 


subject to equation (11) . 
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In order to express Q , and explicitly in terms of P, rewrite 


/\ /\/\/\ 

= F + GPH 
a 


/N \ 

= Ql + H p ' 


/\ 

Ga = + GPI 


(B6) 


where 


F = 


F I 0 
0 I 0 


(Ns+M)x(Ns+M) 


G = 


-U ' 2 

0 I I 


•N 


✓V. 

H 

1 0 


✓N. 

Gw 1 0 

H = 

0 

1 I_ 

(No+M) X (Ng+M) 


0 1 0 . 


(Ng+M) X (Nc+M) 


(Ns+M)x(Nv^+No) 


Qn = 


T I 

Q 

' 2 

0 I 0 


I = 


2 ' h 

0 I 0 


(Ng+M) X (Ng+M) 


-J (No+M)x(n^+No) 


22 = 


-2 ' 2 
0 I 0 


(Nc+M) X (Nc+M) 




(B7) 


One can now write equation (B5) as 


3 j 


= tr 


9p. 


ID 


fK. i 


8P. . 

ID 


Xa + 2 


_ /A. 

9 (GP) 

8p. . 

ID 


(Sx^a) 


+ 2 


3(GP) 

3p. . 

ID 


(iRaGa^A) 
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dJ 









(B8) 


Collecting the derivatives with respect to all the 



results in 


3j 

8p 



. 3-(ax,S-) 



(B9) 


Note that in equation (B9) , both the left and right terms represent 
(Nc + M) X (No + M) matrices. Hence, the evaluation of the trace for 

derivatives with respect to each Pij is no longer required. By expand- 
ing equation (B9) by using equations (12), (19), and (B7) , equations (16) 
to (18) are obtained. 
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SIGNIFICANCE OF A, B, AND C FOR FULL-ORDER CONTROLLER 
For a full-order controller (i.e., M = Ng and R = I) , if we set 


A = F - BH + G^C 


(Cl) 


and use B and C as design variables, then equations (26) can be expressed as 


II = Ac(biv - PH^) + (As/ + Ac)(xg - Xsc)h^ 


(C2) 


II = (q2C + gJs)x^ + g/(a3 + A3^)(x 3^ - x^) 


(C3) 


where 


P I Xg + ^sc ^sc 




S ( Ag + + ^sc 


X) 


(C4) 


(C5) 




e = xc - xs 


(C6) 


Here P is recognized as the steady-state covariance matrix of the estimation 
error e. From equations (11) and (20 ) , one can show that P and S satisfy 
a set of Lyapunov equations as follows. 

Partition equation (11) into four Ng x Ng blocks, i.e.. 




GuC 

A 


Xs_ I X 
I x^ 

'^SC I 


sc 


Xg I X 

;--x . x; 

^SC I 


sc 


BH 


-|T 


GuC 

A 


Gw^w^ 
0 


—WWW I w 


BR^B 


= 0 (C7) 
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Equation (C7) represents four sets of matrix equations: 


FXg 

+ GuOXsc"^ + 

T 

XgF + 

T T 

XgcC Gu 


= 0 

(C8) 


+ G^CXc + 

T T 

XgH B + 

T 


= 0 

(C9) 

BHXg 

T 

T T 

Xsc F + 

T T 
XcC G^ 


= 0 

(CIO) 

BHXgc 

+ AXc + 

T T T 

Xsc H B + 

T 

XcA 

T 

+ BR^B 

= 0 

(Cll) 


Similarly, equation (20) can be expressed in a partitioned form as 


F J 

T 

^s • ^sc 

_ 1 _ 

^S I ^sc 


F_ 1 fuf 

1 

Hd'^QiHd I 0 

BH 1 ~A~ 


“sc 1 

1 

“sc 1 


BH 1 ~A 


0 T 

! ^ ^2^ 


(C12) 


Because of symmetry, equations (C9) and (CIO) represent the same set of equa- 
tions. By substituting equation (Cl) for A and siibtracting both (C9) and (CIO) 
from the sum of equations (C8) and (Cll) , all G^C terms cancel and the follow- 
ing equation is obtained: 


(F - BH)P + P(F - BH)"^ + 


+ br^b"^ 


0 


(C13) 


By substituting equation (Cl) for A and adding all four sets of matrix equa- 
tions represented by equation (C12) , all BH terms cancel and the following 
equation is obtained: 


^F + GuC^'^S + s(f + GuC) + + C^Q2C = 0 (C14) 

From equations (C2) and (C3) , it is easily seen that dJ/dB = 0 and dJ/dC = 0 
provided that 


B 


PH R, 


-1 


Bo 


(C15) 
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C - -Q2 S - Cq 


(C16) 


and provided that 


and 


(^C - ^sc) = 0 
-^sc ) ~ ® 


(C17) 


(CIS) 


which is indeed true as shown below. 

Subtracting equation (Cll) from equation (C9 ) , one obtains 

(F - BH)^Xsc - Xc) + (Xsc - Xcj(p + GuC)'^ + (ph'^ - = 0 (C19) 

Adding the last two sets of partitioned equations (C12 ) , one obtains 

(F - BH)'^^Asc'^ + Ac) + (Asc'^ + Ac)(f + G^c) + c'^^Gu'^S + Q2c) = 0 (C20) 


If equations (CIS) and (C16) are satisfied, then equations (C19) and (C20) are 
homogeneous equations in (^sc*^ ■*" ^c) f respectively, and can have 

only a trivial solution provided that (F -- BH) and -(Ft G^C) have no eigen- 
values in common. Hence equations (C17) and (CIS) are satisfied. Equa- 
tions (CIS) and (C16) for B and C are recognized as the Kalman estimator 
gain and optimal full-state feedback gain matrices, respectively, as given by 
the LQG solution. 

It is also recognized that equation (C17) implies that 


ijeKc"^ = 0 


(C21) 


i.e., the estimation error and the estimated states are uncorrelated, an impor- 
tant result in optimal filtering and prediction. Substituting equations (CIS) 
and (C16) into equations (C13) and (C14) , respectively, results in the steady- 
state Riccati equations for computing P and S : 


FP + PF*^ - Ph'^Rv^^HP + = 0 


(C22) 
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f'^S + SF - SG^Q2 


-1 


G„ S 


+ Hd^QiHd 


0 


(C23) 


Also from equations (C4) and (C17) , the state covariance matrix is given by 


Xg = P + Xc 


(C24) 


and the performance index given by equation (15) becomes 


J = tr 


+ c'^Q2c)xc + ^Hd^QiHd)p 


(C25) 
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DISCUSSION OF ESTIMATION ERROR 


For M = Ng, the states of the controller discussed in this paper are 

optimal estimates of the plant states Xg as shown in appendix C. By analogy, 
for M < Ng, the controller states x^ may be treated as an estimate of M key 
plant states denoted by Xj^. Let the remaining Ng - M plant states be denoted 
by X 2 - Rearranging and partitioning Xg and equations (4) and (5) accord- 
ingly, one can rewrite them as 


X 




(Dl) 


y = 



+ V 


(D2) 


where 


Xi E RX3 


F^^ E RFR^ 


Gui - RG^ 


G = RG 
wl w 


E HR-^ 


Now defining the error e in the partial state estimation as 


e = x^ - x, = x^ - Rx^ 
c 1 c s 


(D3) 


and subtracting the first M equations in equation (Dl) from equation (7) , one 
obtains the error dynamics equation as 


e — ^^11 BH 2 ^)x 2 ^ ^^12 ~ ^^2 ^^2 ~ ^ul^ ” ^wl^ 


(D4) 
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If we let 


A = 


^11 - BHl + 




(D5) 


then equation (D4) becomes 


e = - BH 3 ^)e - (Fj^2 “ ^H 2 )X 2 ~ “ Cx^) - G^jyj + Bv (D6) 

In the usual case^ (u - Cx^) is set to zero as per equation (8) . In the present 
design f to improve the robustness of the reduced-order controller, (u - Cx^) 
is replaced by a fictitious input white noise process of intensity 

(ref. 10) . Thus, 

4 = - (Fj^2 ■ 3^2 )X 2 - Guin^ - + Bv (D7) 


Thus f|u appears as an additional noise parameter in the error dynamics. The 
underlined terms contribute to additional estimation error apart from the white 
noise inputs. These terms are referred to as observation spillover terms in 
reference 16, in which use of a prefilter is suggested to eliminate them. 

If “ ®^1^ stable in equation (D7) and F^ (eq. (10)) is stable 

for a specific choice of B and C, then e is stable and asymptotically con- 
verges to a steady-state value. 
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SYMBOLS 


A 


ac 
. * 


controller dynamics matrix 
dynamics matrix of actuator 
dynamics matrix of turbulence 


Pn] - pm] 

pn]c'[%], 

PJp’PJ, 

Pnl’Pm], 


matrices for polynomial fit (eq. (2)) 
control surface component of Pn] and 
modal component of pn] and pm] 


B 


B 


ac 


ac 

* 


p=] 

e 

F 

Fn.F 


11 '" 12 " • 
* 


gust component of Pn] and pm] 
controller input matrix 

control distribution matrix of actuator 

control distribution matrix of turbulence 

Kalman estimator gain matrix 

controller output matrix 

output matrix of actuator 

output matrix of turbulence 

optimal full-state feedback gain matrix 

reference chord length 

structural damping matrix 

estimation error vector 

plant dynamics matrix 

augmented dynamics matrix 

. components of partitioned matrix F 

dynamics matrix for equation (1) 

matrix defined in equations (B7) 


augmented noise input matrix 


<H 


(eq. (DD) 


^Ul'‘^u2 

^wl'^w2 

/s 

g,g„ 


plant input matrix 
components of partitioned matrix 
gust input matrix 
components of partitioned matrix (eq- (Dl)) 

matrices defined in equations (B7) 


G 


* 

1 




G 


•k 

2 


g 

H 


H 


D 



/N 

H 


* * 
I 


J 

[K] 

L 

I 

M 

[M] 

Nc 

Nd 


N. 


N, 


w 


input matrices in equations (3a) and (A7) 

gravitational acceleration 
sensor output matrix 
design output matrix 

components of partitioned matrix H (eq. (D2) ) 
matrix defined in equations (B7) 

measurement matrices in equations (3b) and (A9) 
unity matrix 

matrix defined in equations (B7) 
performance index 
generalized stiffness matrix 
number of aerodynamic lag terms 
scale of turbulence 
order of controller 
generalized mass matrix 
number of control inputs 
number of design outputs 
number of outputs or sensors 
momber of plant states 
number of plant noise inputs 



p 

/N 

P 




‘‘a 

21 

22 

[q]^[2]c 

Q1/Q2 

q 

‘^f 

R 

^a 

S 

S 

t 

U 

u 

u* 

V 

V 

w 


X ,X ,X 
s' c sc 


Kalman estimator error covariance matrix 

matrix defined in equation (B2) 

nonzero elements of P matrix 

augmented symmetric weighting matrix 

design output symmetric weighting matrix 

control input symmetric weighting matrix 

S“plane approximation of unsteady aerodynamics matrices 

matrices defined in equations (B7) 

dynamic pressure 

flutter dynamic pressure 

key state selection matrix 

augmented noise intensity matrix 

fictitious input-noise intensity matrix 

measurement noise intensity matrix 

plant noise intensity matrix 

full-state feedback optimal solution matrix 

Laplace variable 

time 

input covariance matrix (steady state) 

control input vector 

input vector for equation (1) 

free-stream velocity 

measurement noise vector 

plant noise vector 

covariance matrix of augmented state vector (steady state) 
components of partitioned matrix X^ (eq, (12)) 
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6 

6 


'u 


A 


s c 




- . ,x 


.A 


sc 


augmented state vector 

controller state vector 

augmented plant state vector 

vector of key plant states to be estimated 

vector of plant states not estimated 

state vector for equation (1) 

actuator state vector 

generalized coordinate vector 

time derivative of generalized coordinate vector 
turbulence state vector 
aerodynamic states 
measurement output vector 
design output vector 
accelerometer output 
aerodynamic lag 
control surface deflection 
control surface deflection rate 
commanded control surface deflection 
augmented noise vector 
fictitious input-noise vector 
Lagrange multiplier matrix 

components of partitioned matrix A (eq. (19)) 
generalized coordinate vector for control surface deflection 
generalized coordinate vector for flexible modes 


gust velocity vector 


root-mean-square gust velocity 
wg 

$ matrix of modal amplitudes at sensor location 

Subscript: 

rms root mean square 

Mathematical notation: 

E[ ] expected value of 

tr[ ] trace of matrix 

= defined as 

[ ] matrix 

{ } vector 

T 

[ ] vector or matrix transpose 

[ inverse of square matrix 

Dots over symbols denote differentiation with respect to time. 
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TABLE I.- SUMMARY OF RESULTS FOR 25TH-ORDER PLANT 


|Qj^ = 0.0001, Q 2 = 50000, H = Hp, R^ = 1.0, R^ = 1.0, q = 7.66 kP^ 


Steady- state 


Order Design Optimal responses with Gain margin Phase margin Flutter 

of input noise performance = 0 dynamic 


controller, intensity, 

index. 

• • 

1 

i 1 pressure. 

M Ru 

1 L 

J 

^rms' ^rms^ ' ^rms^ 1 '^rmsf 
deg ■ deg deg/sec ‘g units 

dB 1 rad/sec 

dB : rad/sec deg rad/ sec deg rad/sec 













TABLE II. ^ COMPARISON OF RESULTS WITH CONTROLLERS DESIGNED BY VARIOUS METHODS FOR 25TH-ORDER PLANT 


= 0.0001, Q 2 = 50000, H Hj^, 


R^ = 1.0, R^=1.0, q=7.66 kP^ 


Controller 

design 

method 

Order 

of 

Design 
input noise 

Optimal 

performance 

Steady- state 
responses with 
Ru = 0 


Gain margin 



Phase margin 

controller, 

M 

intensity, 

Ru 

index , 
J 

^rms ' 
deg 

*^rms t 
deg/sec 

^rms ' 
g units 

dB 

rad/ sec 

dB 

rad/sec 

deg 

rad/sec 

deg 

rad/sec 

LQG 

Full 

0.00001 

51.3 

18.0 

795.3 

3.12 

-5.7 

60 

18.9 

616 

-47 

46 

65 

82 

Truncation 

4 

.00001 

110.9 

32.6 

1263.7 

2.39 

-10.3 

59 

6.5 

254 

-17 

38 

57 

129 

Residualization 

4 

.00001 

85.5 

17.1 

665.4 

2.82 

-7.1 

62 

6.0 

164 

-55 

46 

38 

97 

Reoptimi zation 
(present 
method) 

4 

.00001 

59.5 

17.0 

756.0 

3.52 

-5.0 

60 

12.3 

314 

-46 

47 

^53 

83 


























Phase margin 



Point Phase margin 

A -47° at 46 rad/sec 

B +65° at 82 rad/sec 


Gain margin 

C -5.7 dB at 60 rad/sec 
D +18.9 dB at 616 rad/sec 



Figure 5.- Nyquist diagrams of plant plus full-order control law 
(arrows indicate increasing frequency), q = 7.66 kPa. 





ro lO vo 



Real part, rad/sec 

Figure 6.- Closed-loop dynamic-pressure root locus using full-order control law 

= 0 . 00001 . 



Performance index ratio, J/J. 


Cn 

O 



Figure 7.- Typical convergence pattern of performance index using a conjugate 

gradient algorithm. = 0.00001. 







Point Phase margin 2 ” 

A -41° at 52 rad/sec 

B +22° at 92 rad/sec 

Gain margin 

C -4.1 dB at 63 rad/sec 
D +2.6 dB at 120 rad/sec 



(a) = 0. 


Point Phase margin 

A -46° at 47 rad/sec 2 *“ 

B +53° at 83 rad/sec 


Gain margin 

C -5.0 dB at 60 rad/sec 
D +12.3 dB at 314 rad/sec 



(b) = 0.00001. 


Figure 8.- Nyquist diagrams of plant plus fourth-order control law 
(arrows indicate increasing frequency) . q = 7.66 kPa. 


un 
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Dynamic pressure, kPa 

Figure 12.- Variation of steady-state root-mean- square (rms) responses with dynamic pressure 

for fourth-order control law. = 0.00001. 













Phase margin 
-47^ at 46 rad/sec 
+65^ at 82 rad/sec 
Gain margin 
-5.7 dB at 60 rad/sec 


+18.9 dB at 616 rad/sec 
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